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Abstract. In a few years astrometry with the venerable combination of Hubble Space Telescope 
and Fine Guidance Sensor will be replaced by SIM, GAIA, and long-baseline interferometry. 
Until then we remain a resource of choice for sub-millisecond of arc precision optical astrometry. 
As examples we discuss 1) the uses which can be made of our parallaxes of galactic Cepheids, and 
2) the determination of perturbation orbital elements for several exoplanet host stars, yielding 
true companion masses. 
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1. Introduction 

In the early 1970's, Bob O'Dell, then director of Yerkes Observatory, twisted the arm of 
Bill van Altcna, encouraging Bill to propose to study the feasibility of doing astrometry 
with what was then called the Space Telescope (ST). Bill was joined by Larry Frederick 
and Otto Franz. Their efforts insured that astrometry with a Fine Guidance Sensor (FGS) 
would be an integral science component for ST. In 1977 Bill Jcffcrys formed a team to 
respond to the ST Announcement of Opportunity. The first author was a member of 
that team, which then included Raynor Duncombe, Paul Hcmcnway, and Pete Shelus. 
We were eventually chosen as a component (along with van Altena, Franz, and Fredrick) 
of the final STAT (Space Telescope Astrometry Team). Barbara McArthur joined us in 
1985. 

ST was launched in 1990, becoming the Hubble Space Telescope (HST). Because of 
seismic disturbances caused by day/night cycle temperature effects on the original solar 
arrays and slow, but fortunately steady, changes in the FGS, we were unable to prove that 
precision astrometry could be carried out with the FGS until after the first refurbishment 
mission in late 1992. With new, re-designed solar arrays, HST and the two FGS units we 
have used over the years, originally predicted to produce single-measurement astrometry 
with 3 millisecond of arc (mas) precision, have both routinely produced 1 mas or better 
per-observation precision astrometry. Calibration issues, many of which are non-trivial, 
are discussed in McArthur et al. (2002, 2006), Benedict et al. (1994, 2002a), and for 
binary star fringe morphology analysis, Franz et al. (1998). 

The many results from the original STAT include parallaxes of astrophysically inter- 
esting stars (Benedict et al. 2000, 2002a, 2002b, 2003, McArthur et al. 1999, 2001), a 
parallax for the Hyades (van Altena et al. 1997), a link between quasars and the HIP- 
PARCOS reference frame (Hemenway et al. 1997), determination of low-mass binary 
star masses (Franz et al. 1998, Benedict et al. 2001), and searches for Jupiter-mass 
companions to Proxima Cen and Barnard's Star (Benedict et al. 1999). 
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Since 1999, Barbara and I have carried on as General Observers, contributing to the 
study of the lower main sequence mass- luminosity relationship (Henry et al. 1999), the 
intercomparison of dwarf novae (Harrison et al. 2004), and parallaxes of cataclysmic 
variables (Beuermann et al. 2003, 2004; Roelofs et al. 2007) and the Pleiades (Soderblom 
et al. 2005). Jacob Bean, now a postdoc, joined us as a graduate student in 2002. 

Two major themes of our investigations over the last few years have included the 
cosmic distance scale and extrasolar planetary systems. The remainder of this paper 
will outline recent results, specifically drawing from our work (Benedict et al. 2007) 
on the galactic Cepheid Period-Luminosity Relationship (PLR), and the determination 
of extrasolar planetary masses (Benedict et al. 2002c, McArthur et al. 2004, Benedict 
et al. 2006, Bean et al. 2007). 



2. The Galactic PLR 

Our goal was to determine trigonometric parallaxes for nearby fundamental mode 
Galactic Cepheid variable stars. Our target selection consisted in choosing the nearest 
Cepheids (using HIPPARCOS parallaxes, Perryman et al. 1997), covering as wide a 
period range as possible. These stars were the brightest known Cepheids at their respec- 
tive periods. Our new parallaxes provided distances and ultimately absolute magnitudes, 
M, in several bandpasses. Additionally, our investigation of the astrometric reference 
stars provided an independent estimation of the line of sight extinction to each of these 
stars, a contributor to the uncertainty in the absolute magnitudes of our prime targets. 
These Cepheids, all with near solar metallicity, should be immune to variations in ab- 



solute magnitude due to metallicity variations, e.g. Groenewegen et al. (2004) Maori 



et al. (2006) Adding our previously determined absolute magnitude for 5 Cep, Benedict 
et al. (2002b) | we established V, I, K, and Wyi Period-Luminosity Relationships using 



ten Galactic Cepheids with average metallicity, ([Fe/H])=0.02, a calibration that can be 
directly applied to external galaxies whose Cepheids exhibit solar metallicity. Details are 



given in Benedict et al. (2007) 



One of our interesting results is an independent determination of the distance mod- 
ulus of the Large Magellanic Cloud (LMC) , particularly given the recent exploration of 
possible socialogical influences on such determinations (Schaefer, 2007). OGLE (Udalski 
et al. 1999) has produced the largest amount of LMC Cepheid photometry. We worked 
with an apparent Wyi PLR for 581 Cepheids in the LMC. These data were carefully 
preened, selecting only Cepheids with normal light curves and amplitudes from |Ngeow| 
et al. (2005) They provide a zero-point of 12.65 ± 0.01. Direct comparison of that Wyi 



zero-point with the zero-point of the left-hand M H -(\v) PLR. in Figure j 1 1 yields an LMC 
distance modulus 18.51 ±0.03 with no metallicity corrections. 



Macri et al. (2006) demonstrate that a metallicity correction is necessary by comparing 



metal-rich Cepheids with metal-poor Cepheids in NGC 4258. With a previously measured 
[O /H] metallicity gradient from Zaritsky et al. (1994) [ Macri et al. find a Cepheid metal- 
licity correction in Wy;, 7 = —0.29 ± 0.09 r ± 0.05 s magnitude for 1 dex in metallicity, 
where r and s subscripts signify random and systematic. This value is similar to an ear- 
lier Wvi metallicity correction (Kcnnicutt et al. 1998) derived from Cepheids in M101 
(-0.24± 0.16). Taking the weighted mean of the Kcnnicutt and Macri values and using 
the difference in metallicity of LMC and Galactic Cepheids (-0.36 dex from means of 
the data in Groenewegen et al. 2004 tables 3 and 4) we find a metallicity correction of 
-0.10± 0.03 magnitude with the Galactic Cepheids being brighter. With this estimated 
metallicity correction we obtain a corrected LMC modulus of 18.41 ± 0.05. One other 
recent determination is noteworthy for its lack of dependence on any metallicity correc- 
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tions. 

close agreement with our new value. 



Fitzpatrick et al.(2003) derive 18.42± 0.04, from eclipsing binaries, a modulus in 
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Figure 1. The densification of the PLR. Our original Cepheid parallaxes (left) have calibrated 
techniques applied to many other, more distant Cepheids (right, after Fouque et al. 2007). Note 
the reduction in the error of the slope of the PbR. 

As an additional check on our adopted LMC modulus, m-M= 18.41 ± 0.05, a Macri 
et al. (2006) differential modulus (LMC-NGC 4258) leads to a modulus of 29.29 ± 0.08 
for NGC 4258, a value in excellent agreement with the maser-based distance for NGC 
4258, m-M = 29.29 ± 0.15 from |Herrnstein et al. (1999) 



One of the first uses of our new Cepheid parallaxes was the calibration of various tech- 
niques used to provide distances for many more long period (and all far more distant) 



Cepheids as described in |Fouque et al. (2007) These techniques include infrared surface 
brightness and interferometric Badde-Wesselink, which compare measured physical di- 
ameters with radial velocities. This calibration results in a densification and extension 
of the Galactic Cepheid PLR, and is illustrated in Figure [T] The immediate value of 
the densification is a greater confidence in the absolute magnitudes of the longer-period, 
brighter Cepheids, the ones critical for extragalactic distance determination. 



3. Extrasolar Planet Masses 

Currently, fewer than 10% of the more than 200 candidate exoplanets 



(http://exoplanet.eu) orbiting nearby stars have precisely determined masses. Because 
the most successful technique for detecting candidate exoplanets, the radial velocity 
method, suffers from a degeneracy between the mass and orbital inclination for most of 
the known exoplanet candidates, only a minimum mass, Msini, is known. In principle, 
radial velocities alone can be used to determine the masses of exoplanets in multiplanet 
systems, if two or more planets are experiencing significant mutual gravitational interac- 
tions on short timescales. However, only one such system has been investigated, results 
from different groups vary significantly, and the effects of non-coplanarity remain to be 
considered. 

Thus, establishing precise masses, rather than arguing statistically (using Msini, which 
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everyone does, e.g. the otherwise excellent review of Udry & Santos 2007), for the majority 
of exoplanet candidates requires additional observations. The primary techniques that 
have been employed to break the mass-inclination degeneracy in radial velocity data are 
astrometry (Benedict et al. 2002c, McArthur et al. 2004, Benedict et al. 2006, Bean 
et al. 2007) and transit (e.g., Bakos et al. 2007) observations. Ideally, all candidates 
should have their masses determined. Obviously, not all systems are edge-on to us. And, 
unfortunately, the size of the perturbation due to the companion decreases linearly with 
distance, placing many of them beyond the reach of HST FGS astrometry. 



Table 1. Masses from HST FGS Astrometry 



Component 


M,[M ] 


[Fe/H] 


d(pc) 


ecc. 


M p {M Jup ) 


a (mas) 


inc. (°) 


P (d)| 


GJ 876 b 1 


0.32 


-0.12 2 


4.7 


0.1 


1.9±0.5 


0.25 


84±6 


61 


55 Cnc d 3 


1.21 


0.32 


12.5 


0.33 


4.9±1.1 


1.9 


53±7 


4517 


e Eri b 4 


0.83 


-0.03 


3.2 


0.7 


1.6±0.2 


1.9 


30±4 


2502 


HD 33636 B 5 


1.02 


-0.13 


28.1 


0.48 


142±11 


14.2 


4.1±0.1 


2117 



Benedict et al. 2002c, 2 Bean et al. 2006a, 2006b, 3 McArthur et a/.2004, 4 Benedict et al. 2006, 
5 Bean et al. 2007 

To partially rectify this lack of knowledge of the true masses, the HST allocation 
process awarded us over 100 orbits in Cycle 14 and 15 to be used to determine actual 
masses for HD 47536 b, HD 136118 b, HD 168443 c, HD 145675 b, HD 38529 c, and 
HD 33636 b. The object chosen for discussion in this contribution, HD 33636 (Bean 
et al. 2007) presents a cautionary tale. As we shall see, at least in some cases, Msini is 
not a very good estimator of the actual companion mass. 

We chose to work first on HD 33636 because we had finished collecting all the HST data 
for this object. Radial velocity data (Figure 2) were obtained from past investigators (Lick 
and Keck data from Butler et al. 2006 and Marcy (2007, private communication) , Elodie 
data from Perrier et al. (2003)) and from a high-cadence campaign with the Hobby- 
Eberly Telescope High-Resolution Spectrograph (HRS; Tull 1998). The latter data were 
used to search for shorter-period companions, the unknown and unmodelcd presence of 
which would have introduced noise into our velocity data. No such companions were 
found with detection limits detailed in Bean et al. (2007). Limits depend on period 
and assumed eccentricity of the tertiary. For example these velocity data would have 
permitted detection of a tertiary with M > 0.5Mj up for P < 300 d and eccentricity 
< 0.7. 

Our techniques for combining radial velocity and HST FGS astrometry have been 
detailed in Benedict et al. (2001, 2006) and Bean et al. (2007). For HD 33636 we found 
an inclination, i = 4.0° ± 0.1° and, as shown in Figure 3, a perturbation semimajor axis, 
a = 14.2 ± 0.2 mas. Assuming a mass for the primary, HD 33636, M = 1.02 ± 0.03 
M Q (Takeda et al. 2007), we obtain a mass for HD 33636 b, M = 142 ± llM Jup = 
0.14 ± 0.01 M Q . HD 33636 b is actually HD 33636 B, an M6V star. Including it in any 
statistical studies depending on Msini, such as metallicity, eccentricity, and semimajor 
axis, would corrupt those relationships. 

We have made some progress. Table 1 collects the mass estimates that have already 
resulted from our combination of HST astrometry and high-precision radial velocities. 
Over the next year we will finish our analysis and determine actual masses for HD 47536 
b, HD 136118 b, HD 168443 c, HD 145675 b, and HD 38529 c. But, before those results, 
we will have finished our analysis of the v Andromedae system (McArthur et al. 2008) . 
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We have a preliminary indication that for the first time the degree of coplanarity of 
an extrasolar planetary system associated with a normal, main sequence star will be 
established. Beyond that, we are in the process of acquiring HST FGS data for similar 
coplanarity tests on the multiplanet systems associated with HD 128311, HD 202206, 
/i Ara, and 7 Cep, paving the way for future studies of hundreds of such systems with 
ground-based long-baseline interferometry and the space-based projects, SIM and GAIA. 
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Figure 2. HD 33636: Radial velocities and residuals to our final orbit from combined 
velocities and astrometry. The HET residuals are slightly better than the Keck residuals. 
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Figure 3. The astrometric perturbation due to the stellar companion to HD 33636. Based 
on radial velocities only this object had a minimum mass Msini — 9AMj up . The inclination, 
i = 4.1°, and size of the semi-major axis of the perturbation, a = 14.2 mas, identify the 
companion as an M dwarf star. The box marks the date of the next periastron passage. 
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